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Series and e-expansion of the hypergeometric functions 
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Recent progress in analytical calculation of the multiple {inverse, binomial, harmonic} sums , related with 
e-expansion of the hypergeometric function of one variable are discussed. 
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1 . In the framework of the dimensional regular- 
ization [1] many Feynman diagrams can be writ- 
ten as hypergeometric series of several variables 
[2] (some of them can be equal to the rational 
numbers) . This result can be deduced via Mellin- 
Barnes technique [3] or as solution of the differ- 
ential equation for Feynman amplitude [4] . How- 
ever, for application to the calculation of radia- 
tive corrections mainly the construction of the e- 
expansion is interesting. At the present moment, 
several algorithms for the e-expansion of different 
hypergeometric functions have been elaborated. 
They were mainly related to calculations of con- 
crete Feynman diagrams [5]. Only recently, the 
general algorithm for integer values of parameters 
has been described in [6] and its generalization 
has been done in [7]. The results of expansion 
are expressible in terms of the new functions, like 
harmonic polylogarithms [8] or their recent gener- 
alization [6,9]. Let us shortly describe, how this 
algorithm woks on the example of the general- 
ized hypergeometric function of one variable. The 
starting point is the series representation: 
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where (a)j = T(a + j)/T(a) is the Pochham- 
mer symbol. We concentrate on the case q + \Fq, 
when series converges for all \z\ < 1, and on the 
integer or half-integer values of the parameters 
{Ai, Bj} 6 {rrii,mj + 5}. To perform the e- 
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expansion we use the well-known representation 

(to + ae)j 
(m)j 

exp|-f;^^[5 fe (m+i-l)-5 fc (m-l)]| , 

where to is an integer positive number, to > 1 and 
Sk(j) — Ya=i l~ h 1S the harmonic sum satisfying 
the relation S k {j) = S k (j - 1) + l/j k . For half- 
integer positive values Ai = rrii + 1/2 > 0, we use 
the duplication formula 
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To work only with positive values for parameters 
of hypergeometric function we can apply several 
times the Kummer relation: 
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After applying this procedure the original hyper- 
geometric function can be written as 
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with 
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In this way, the e-expansion of the hypergeomet- 
ric function (1) is reduced to the calculation of 
the multiple sums 
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where {rrij, n k , l w ,Pr} - positive integer numbers. 
In the calculation of massive Feynman diagrams 
[10-12] we get multiple sums of the following 
form, 
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where u is an arbitrary argument and we accept 
that the notation S a and Sb will always mean 
S a (j — 1) and Sb (2 j — 1), respectively, even we do 
not mention this explicitly When there are no 
sums of the type S a or §b in the r.h.s. of Eq. (2), 
we put a "— " sign instead of the indices (a) or (b) 
of S, respectively. Some indices (a) or (b) may be 
equal to each other, which is equivalent to power 
of a proper harmonic sum. For particular values 
of k, the sums (2) are called 



harmonic 

k = < 1 inverse binomial 

-1 binomial 



sums 



These sums are related to e-expansion of the hy- 
pcrgeometric functions of type (1) with the fol- 
lowing set of parameters: 
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u = A k z . (3) 



In the recent paper [12], the sums of type (2) up 

to weight 4 have been studied in detail. 

2. Let us rewrite the multiple sums (2) in the 

following form, Y> {k) B . c {u) = Y^Li uJ Va;B;c(J) > 
where A = (oi,...,a p ) and £? = (&i,...,6 g ) 
denote the collective sets of indices, whereas 
Va-b-cU) i s tne coefficient of v? 

v'aUj) 
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The idea is to find a recurrence relation with re- 
spect to j, for the coefficients Va-b-cU) ano - then 
transform it into a differential equation for the 
generating function T, A B . c (u). In this way, the 
problem of summing the series would be reduced 
to solving a differential equation [13]. Using the 

explicit form of t)a-b- c (J) gi ven m Eq. (4), the re- 
currence relation can be written in the following 
form: 

[ 2 (2j+i)] fc (j+irsi% ;c (j+i) 

= 3 C VA k Lc(3) + r%(j), (5) 

where the explicit form of the "remainder" 

r i;s(i) is s iven b y 
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In other words, it contains all contributions gener- 
ated by j~ afc , (2j)~ 6; and (2j+l)~ b < which appear 
because of the shift of the index j. Multiplying 
both sides of Eq. (5) by u J , summing from 1 to 
infinity, and using the fact that any extra power 
of j corresponds to the derivative u(d/du), we ar- 
rive at the following differential equation for the 
generating function T, A . B . c (u): 
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where i2^/ B (w) = X)^=i wJ 'm-b (■?')■ The r.h.s. of 
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differential equation for sums includes the multi- 
ple sums with shifted index 
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where we accept the same conditions for the in- 
dices {a{\ and {bj}, as in the previous case. For 
investigation of these sums we again apply the 
generating function approach. In this case, the 



recurrence relations for the coefficients v 
are the following: 
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with r^ B (j) given by Eq. (6). The proper set of 
differential equations is 

= S P o + {u£ + l)R%(u), (12) 
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Equations, [(7), (12)], [(8), (13)], [(9), (14)], form 
the closed system of differential equations. 
3. From the analysis, given in [12] we have de- 
duced that the set of equations for generating 
functions has a simpler form in terms of a new 
variable. For multiple inverse binomial sums it is 
defined as 
\/u~- 
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and for multiple binomial sums it has the follow- 
ing form: 

1 -y/l-Au X 



X = , u = — -7T. (16) 

l + vT=4ir (l + x) 2 

Let us consider the differential equation for mul- 
tiple inverse binomial sums in terms of new vari- 
ables. Equation (7) takes the form 
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Equation (17) could be rewritten as 
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or, in equivalent form: 
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From this representation we immediately get the 
following statement: 

If for some j the series £ .. B -(y) are express- 
ible in terms of harmonic polylogarithms, the 

sums ^A-B-j+i(y) can a ^ so be presented in terms 
of harmonic polylogarithms. This follows from 
the definition of the harmonic polylogarithms (see 
Ref. [8]). 

In a similar manner, let us rewrite the equation 
for generating function of the multiple binomial 
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which could also be rewritten as 
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or, in an equivalent form: 
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Again we get the previous statement. 
4. The differential equation for multiple inverse 
binomial sums with the shifted index has a more 
complicated form. For their analysis let us use the 
geometrical variable [14] defined via ug = 4 sin 2 | 



(0 < ug < 4). In terms of this variable, Eq. (12) 
could be written as 
l 
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This equation can be decomposed into the system 
of differential equations 
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The formal solution of Eq. (28) is 

gA;B;c(0) = \ Sin 6R% B ( U0 ) 
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where u^ = 4 sin ^ . Substituting this result in 
Eq. (27) and integrating one time we get 
d 
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(30) 



For c = 1 the r.h.s. of Eq. (30) divided by sin | is 

the solution for the sum G A : B . 1 (ug). Let us now 

apply this approach to the sum G_._. 3 (ug) which 
was not solved explicitly in Rcf. [12]. Using the 
relation (see Eqs. (27)) 



3-fe 

2tan|— ) PA;B;c(0) = PA:B:k{0) 
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and expression for /9_ ; _ ;2 (#) (see Eq. (2.74) in 
[12]), p-;_ ;2 (0) = 2Ti 2 (tan |) - sinf , we get 
after integration by parts 
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smfj-sinf-f ^-2^J 



= k[p_;-; 2 (0) + 
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where we have used the short-hand notation 
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the integral 1(9) is defined as 
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is the generalized log-sine function [15]. In partic- 
ular, Ls^ (6) = Lsj(9). The integral (33) can be 
evaluated in terms of the polylogarithms of the 
complex argument [15] 

1(6) = -§< 2 -3Ls 3 (^ - 8)-±Ls 3 (29)+Ls 3 (9) 
+41n2[Ls 2 (7r-60 + Ls 2 (l9)] + 83?Si, 2 (itan§) 
+29 In (tan f) [in 2 -In (cos f)] , (35) 

where 
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sm £ — 



1 /)3 



8^ iz V""' 2/ 48" ' (^^) 

and Ti 3 (x, a) is the inverse tangent integral of 
two variables. Collecting all expressions we get 

P-.,-., 3 {0) = §Li 2 (sin 2 f)-4Ti 3 (tan f, cot f) 

+Ls 3 (7r-f)-Ls 3 (§)+iLs 3 (0)-Ls 3 (7r) 

+ln(tanf)[Ls 2 (f)+Ls 2 (7r-f)]+g-sinf 
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Combining this result with Eqs. (2.61), (2.62) and 
(2.54) from Ref. [12] we write the one-fold integral 
representation for the sum S._. 3 . 1 ('U0) 

£-;3;lM = g tan f 
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where Ls 3 (9) and LS4 (9) are expressible in 
terms of Clausen functions (see Appendix A of 
Rcf. [11].) To obtain a result valid in a different 
region of variable ug (—4 < ug < 0), the analyt- 
ical continuation procedure, firstly described in 
Ref. [16] (see also [11]) should be applied. 



5. Our study [12] allows us to construct some 
terms of the e-expansion of the generalized hyper- 
geometric function p + \Fp and obtain new analyt- 
ical results for higher terms of the e-expansion of 
the two- and three-loop maser-integrals entering 
in different packages [17,18]. 

It was noted, that individual terms of our con- 
struction include the log-sine functions of argu- 
ments 6</2 and it - 9/2 (see Eqs. (37, 38) and 
Eq. (2.90) in Rcf. [12]) which, however, are can- 
celled in the considered order of e-expansion of 
the Feynman diagrams. For single scale massive 
diagrams, when 6 = 7r/3, the arguments of func- 
tions are equal to 7r/6 and 57r/6, respectively. It 
opens the question about possible generalization 
of the "sixth root of unity" basis, introduced in 
Ref. [19] (for the recent progress see [11,20]). 
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